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\ We calculate the volume dependence of the ground-state energy of n identical bosons with short- 

£SJ . range repulsive interactions in a periodic spatial volume of side L, up to and including terms of 

order L~ 6 . With this result, Lattice QCD calculations of the ground-state energies of three or more 
pions will allow for a systematic extraction of the three-pion interaction at this order in the volume 
expansion. 
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I. INTRODUCTION 



A major goal of the field of strong-interaction physics is to determine the spectrum of hadrons and nuclei from 
Quantum Chromodynamics (QCD). Lattice QCD is the only known way to rigorously compute strong-interaction 



quantities, and as such, an increasing effort is being put into understanding the lattice QCD calculations that will 
be required to extract even the most basic properties of light nuclei. Ground-state hadron masses, including those 
of nuclei, can be extracted from the long-time exponential behavior of a two-point correlation function with the 
appropriate quantum numbers in Euclidean space. Excited states that are stable against strong decays, such as the 
low-lying energy levels of nuclei, can be extracted from the time- dependence of correlation functions over a large time 
interval with various methods, such as those described in Ref. [H,12|. By contrast, excited states that are unstable 
against strong decays can be extracted by analyzing the volume dependence of the energies of the scattering states 0] . 
However, such direct calculations of the masses of nuclei beyond ^4^6 appear computationally prohibitive using 
t-H ' current lattice techniques. Consequently, the properties of nuclei, beyond the lightest few, will likely be systematically 
f — [ , determined by matching lattice QCD onto nuclear effective field theories (NEFT's) 1 . Lattice QCD calculations of 
few nucleon systems will be used to constrain the low-energy constants (LECs) appearing in the NEFT, which will 
t***'- \ then be used to compute the properties of larger nuclei. To extend calculations to larger nuclei, further matchings (for 
example to the no-core shell model Q) will be required. The most important LECs relate to the interactions in the 
two-body sector, characterized by the effective-range expansion of the two-particle phase shift. These contributions 
(at the physical pion mass) are well constrained by experiment, but can also be determined usin g la ttice QCD. 
Recent work has used the volume dependence [1, [l(| of two nucleon [llj and hyperon-nucleon 12. Il3l | systems to 
extract the corresponding two-body low-energy phase shifts at specific momenta from fully-dynamical mixed-action 
lattice calculations (domain-wall valence quarks on fourth-rooted staggered configurations generated by the M1LC 
collaboration). In principle, phase shifts for all momenta below inelastic thresholds can be determined from QCD 
given sufficient computational resources. 

In order to perform NEFT computations in systems of more than two nucleons, A > 2, it is necessary to include the 
multi-nucleon interactions, such as the three-body force. The leading, local three-body interaction is a momentum- 
independent contact operator, Q 3 - hod y ^ D(fi) [ N N ] 3 , where \i is the renormalization scale (which we will discuss 
subsequently), and N is the nucleon annihilation operator. A determination of the coupling -D(/i) from lattice QCD 
will require the study of three-nucleon systems. This is a complicated task, and aside from the numerical complexities 
of such a lattice calculation, theoretical developments arc required in a number of relevant areas. 

In this work, we compute the ground-state energy of n identical bosons in a finite cubic volume with short-range 
repulsive two-body interactions, subject to periodic boundary conditions. This system was first investigated for hard 



1 For reviews of NEFT see Refs. [1|E@|. 
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spheres by Huang and Yang [lj], and we extend their results to include order L~ 6 contributions in the volume 
expansion, where the two-body effective range and three-body interactions first enter. For the case of three hard- 
spheres, our result agrees with that of Wu [15| (modulo renormalization issues and the appearance of a three-body 
interaction) . While our results are not directly relevant to the study of nuclei, they are necessary for extracting the 
three-boson interaction in the study of multi-boson systems (for instance multi-pion systems). 



II. THE GROUND STATE ENERGY OF n BOSONS 



Motivated by interest in performing lattice QCD calculations, we calculate the ground-state energy of n bosons 
of mass M confined to a finite volume of size L 3 with periodic boundary conditions 2 in non-relativistic quantum 
mechanics. In their classic 1957 paper, Huang and Yang [3] (see also Ref. [HI EH) considered this problem formulated 
for hard spheres using the technique of pseudo-potentials. They calculated the dependence of the n-boson ground 
state on the volume up to L~°: 
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where the sums extend over all three-vectors of integers, and [Z)=nl/(n — k)\/k\. In eq. ([T]), the notation has been 
modified, and numerical values of X and J have been corrected, compared to the expressions presented in Refs. 14,[l6l|. 
The leading term in eq. (P) was derived by Bogoliubov [l7|. In the special case of two particles, Liischer 0, |l0| has 
shown that the energy shift determined in non-relativistic quantum mechanics, given by eq. |T]) with n — 2, can also 
be derived from quantum field theory, and thus is a general result 3 . 

At O (L -5 ), the ground-state energy in eq. {1} is not sensitive to three-body interactions. On dimensional grounds, 
three-body contributions should first enter at O (L~ 6 ) and in this work we extend the result given above to that 
order. This will allow for the extraction of three-body interactions from lattice QCD calculations. The ground-state 
energy of the n-boson system is calculated with an interaction of the form 



V( ri ,...,r n ) = V J2S (3) ( 
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where the ellipsis denote higher-body interactions that do not contribute at the order to which we work (in general, 
m-body interactions will enter at 0(L 3 ( 1 ~ m ^). For an s-wave scattering phase shift, 8(p), the two-body contribution 
to the pseudo-potential is given by r\ = — fjP -1 tan<$(p) = jja + j^-a 2 r(p 2 + p 2 ) + . . ., keeping only the contributions 
from the scattering length and effective range, a and r, respectively. At 0(L~ 6 ) the coefficient of the three-body 
potential, 773, is momentum independent. While up to this point the discussion has been phrased in terms of m-body 
pseudo-potentials, the modern language with which to describe these interactions and calculations is that of the 
pionless EFT, EFT(^) [22|, HH, [24| . In EFT(tF) the divergences that occur in loop diagrams can be renormalized order- 
by-order in the expansion, preserving the power counting. This is somewhat less obvious when using the language of 
pseudo-potentials. 

At O (L~ e ) (order V 4 , corresponding to four insertions of the two-body potential given in eq. ([3])) the energy- 
shift has contributions with interactions amongst up to six particles. Using standard techniques of non-degenerate 
perturbation theory (or the perturbative expansion of EFT(/f)), the ground-state energy of n bosons with repulsive 



2 The effects of the finite time extent of lattice calculations are exponentially small in the temporal length and are i gno red herein. 

3 The two-body result has been extended to the situation where the center-of-mass is moving relative to the volume [I8t [2(| , and also 
to the case in which there are coupled channels [2l| . 
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two-body interactions is 
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which can also be written as 
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In the case of n — 3, the logarithmic behavior with volume was first derived by Wu 15J, and the ground-state energy 
has recently been computed to O (L~ 7 ) by S. Tan [25|. The additional integer sums that contribute at 0(L~ 6 ) are 
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The sums Q and 7Z arise from intrinsic two-loop contributions to the energy and are naively divergent (TZ contains 
a nested divergence, requiring the subtraction to preserve the one-loop scattering amplitude). As is the case for Z, 
given in eq. we define (regulate) these sums using dimensional regularization. It is straightforward to show that 
the dimcnsionally-regulated sums that appear in eq. (0) and eq. (|SJ) are 



TZ -> TZ - 2v / 3tt 3 log(/i L) 
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Numerically, we find that Q = -100.75569 and TZ = 19.186903. In eq. © and in what follows from it, the MS-scheme 
is used to define the three-body coefficient, 773 (/j) in which only the pole, ~ jz^, is subtracted from the divergent 
two- loop diagrams (explicitly absorbed by 773 (/j)). This renormalization, and the associated logarithmic dependence 
on fj, L in eq. (HJ, make clear that the three-body interaction, 773 (fi), must enter at this order [26|; the scale dependence 
of 773 (/_i) exactly cancels that of the logarithms in eq. ^ and in subsequent expressions. The logarithms of \x L that 
appear in eq. (|4]) result from the relation between the momentum-space integrals in <i-dimensions and the integer 
sums. It is worth stressing that the quantities Q and TZ are subtraction-scheme dependent. 

In the first form of our main result, eq. ((3]), the combinatoric factors make clear the number of particles involved 
in the various contributions; terms proportional to (") have j interacting particles and n — j spectators. The second, 
simpler form of our result, eq. follows from evaluating the combinatoric factors. As can be seen from these 
expressions, at this order the two-body effective range also contributes. 
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For two, three and four particles the ground-state energies are 
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From eq.{4|, it is possible to construct combinations of these energies that are directly sensitive to the three-body 
contributions at L~ 6 . That is, for n > 2, 
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This can also be accomplished by forming judicious combinations of the energies of three different systems, such as 
Eq(4, L), Eq(3, L) and Eo(2, L). Further, combinations of energies can be formed in which the three-body interaction 
is absent, such as 



E (2,L)-^E (3,L) + Ie (4,L) 
3 6 
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Relations such as that of eq. (fT5|) will provide a useful check of both the statistical and systematic uncertainties 
associated with lattice QCD calculations. 

The calculation of ground-state energies described here has been derived in a non-relativistic framework, however 
the results remain valid relativistically. In the two-body case, this has been shown by Luscher [?|. In the higher-body 
case, the non-relativistic calculation will not correctly recover a field-theoretic calculation, due to relativistic effects in 
multiple, two-body interactions involving three or more particles. At O (L~ 4 ) , only two-particle interactions contribute 
to the n-body ground-state energy and the results of Ref. @ follow without modification. Since the interaction of 
three particles due to the two-body interaction first enters at L~ 5 , and relativistic effects in such interactions are 
suppressed by (ML) -2 , the first relativistic effects will occur at Q(L~ 7 ). Generally, relativistic effects can be included 
perturbatively into the volume expansion by using EFT(f) [22l. l23l. 1241] . 

Beyond 0(L~ 6 ) a number of other effects are also important: 

• Higher partial waves will contribute to the ground-state energy at finite volume. The cubic periodic boundary 
conditions imposed on the system are such that the spatial symmetry group is H(3) with the ground-state 
wavefunction transforming in the representation (see Ref. [27|). Contributions to the ground-state energy 
from interactions other than s-wave (as classified in the continuum) first enter at order (D(L~ 7 ) via three two- 
body interactions, two of which are s-wave and one of which is p-wave. The leading two-body contribution from 
higher partial-waves results from a single insertion of an / = 4 operator, and enters at 0(L~ n ). 

• It is straightforward to show that the leading, momentum-independent m-body operators contribute to the 
ground-state energy at i 3 ( 1_m ), and therefore, the four-body local operator will first enter at order L~ 9 . It is 
obvious that the momentum-dependent m-body operators first enter two orders higher than the momentum- 
independent operators. For instance, contributions from the two-body shape parameters (higher-order terms in 
the expansion of pcot5(p)) first occur at 0(L~ & ). 
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• In the calculations that we have performed, the symmetrizations that are required to satisfy Bose statistics are 
trivial. At higher orders in the expansion, where multiple insertions of the two-body effective-range interaction 
occur, symmetrization of the intermediate states will have non-trivial consequences and calculations will become 
increasingly complicated. 

III. DISCUSSION AND CONCLUSION 

In this work we have computed the volume dependence of the ground-state energy of n bosons in a cubic volume 
with periodic boundary conditions. Knowledge of this dependence is necessary in order to establish a connection 
between lattice QCD calculations of n-boson systems in Euclidean space and the multi-body interactions contributing 
to the properties of many-body systems in Minkowski space. We conclude by highlighting a number of issues and 
possible extensions of this work: 

• Our results are easily generalized to asymmetric volumes. In the asymmetric case, the integer sums that appear 
in the expressions for the ground-state energy will become dependent on the asymmetry parameters (2a . [29j 
(for some asymmetries, one or more of these sums can vanish [2^|). Further, the reduced symmetry of the 
system forces derivative interactions and higher partial-wave interactions to enter at lower orders in the volume 
expansion [30]. 

• Calculation of excited-state energies are straightforward and tedious. Many contributions (diagrams in the 
perturbative expansion) that vanish in the calculation of the ground-state energy will no longer vanish, and the 
combinatoric factors become somewhat more complicated because of the symmetrization requirements. 

• Our calculation of the ground-state energy of the n-boson system is an expansion whose terms depend parametri- 
cally on n and L as n a /L 13 , where a and (3 are positive integers. Much of the discussion in the literature regarding 
collections of bosons concerns Bose-condensed systems, with the results that are summarized in Ref. [3l|. It 
is clear that in the large-n limit, relevant to the Bose-condensed systems, this expansion fails, and as such, no 
direct connection between our calculation and Bose-condensed systems can be made. 

• We have restricted ourselves to the case of repulsive interactions among the bosons. The reason for this is that 
for an attractive two-body interaction of sufficient strength, the ground-state will not be perturbatively close to 
n non-interacting bosons but will, in fact, be a system of two-body (or higher body) bound states, e.g. one can 
imagine that for the three-body system the ground-state is a two-body bound state interacting with the third 
boson. The finite-volume behavior of the two-particle system with shallow (compared to the inverse range of 
the underlying interaction) bound states was discussed in Ref. [32| . and subsequently numerically explored in 
Ref. [33]. In such systems, the leading finite volume corrections have both power-law and exponential volume 
dependences, and are somewhat more complicated to analyze. 

• The expressions we have presented in this work form the basis for investigating n n correlation functions with 
isospin J = I x = n in lattice QCD. By forming well-defined ratios of correlation functions, as given in eq. (|14[) . 
the three-body interaction can be calculated. The natural choice of renormalization scale is the ultra-violet 
cutoff of the low-energy theory, which is approximately the range of the interaction, u ~ r . This gives rise to 
a large logarithm, log(rL) / L e , which should formally dominate the L~ e contribution to the ground-state energy. 

• As the lattice-QCD study of nuclei is the underlying motivation for this work, it is worth considering difficulties 
that will be encountered in generalizing the result described here to fermionic systems. A significant difference 
between the three-neutron system and the three-boson system is that the Pauli-exclusion principle forces the 
unperturbed ground state of three neutrons to have non-zero momentum. One of the neutrons will have mo- 
mentum |p| = 2ir/L, which provides power-law volume dependence to the ground state even in the absence of 
interactions. As a result the momentum-dependent interactions will enter at lower orders than in the boson 
system. The lowest-energy state with zero total momentum consists of one neutron at rest, and the other two 
moving back-to-back, each with |p| = 2ir/L. 
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